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Abstract. We generalize earlier studies on the Laplacian for a bounded open domain fl £ R with connected 
complement and piecewise smooth boundary. We compare it with the quantum mechanical scattering operator for the 
exterior of this same domain. Using single layer and double layer potentials we can prove a number of new relations 
which hold when one chooses independently Dirichlet or Neumann boundary conditions for the interior and exterior 
problem. This relation is provided by a very simple set of ^-functions, which involve the single and double layer 
potentials. We also provide Krein spectral formulas for all the cases considered and give a numerical algorithm to 
compute the ^-function. 
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1. Introduction 

In an earlier paper [EP2], we derived an identity between the integrated density of states for 
the eigenvalues for the Laplacian in a domain O with Dirichlet boundary conditions, and the 
scattering phases for the exterior of the same domain, also with Dirichlet boundary conditions. 
In the present paper, we derive similar identities, and prove several of them, for the case where 
the boundary conditions can also be of Neumann type. We will discuss and illustrate similarities 
and differences between the various cases. 

Although it is not possible to really formulate the identities we are going to derive without 
making precise definitions, we summarize here the main results in an informal way. We 
consider a bounded domain Qui R 2 , and we let T be its boundary. We denote by G z the free 
Green's function G z = (— A — This is an analytic function of z, except for a logarithmic 
singularity at the origin. We assume that the branch cut is along R + . We let G z (x, x 1 ), with 
x, x G R 2 , denote the integral kernel of G z . Then we define 

A z = G z\rxr > B z = ~ G z N -V lr x r > C z = ~ N ' ^G Z N ' v lr x r > 

where the normal TV points out of il and | rxr is the restriction to the boundary. This is a more 
precise notation for the "normal derivative on the boundary." Setting E = —1, we define 4 
(■-functions, for z £ R + , 

C DD (z) = det {A-]A Z ) , 
C DN (^) = det{{\-B Eo )-\\-B z )) , 
C ND (z) = det{{l + B Eo )-\ l 2 +B z )) , 
( NN (z) = det (C^C Z ) . 

The symbols D and N stand for Dirichlet and Neumann boundary conditions. Then we have 
the following nice identities between the integrated densities of states iV B for the — A in the 
interior and the scattering phases B , for the exterior fl c of il, 1 with boundary conditions 
B.B' e {D.N}. 

tvN b (E) = e B ,(E)-ImlogC B) B'(£ + iO). (1.1) 

Furthermore, 

|ImlogC B ,B'(-E + *0)| < const. E 1/2 log E . (1.2) 
We also have a Krein trace formula, valid for any nice function F: 

Tr(F(-A QjB ©-A QC!B ,)-F(-A)) 

n 

1 We define the scattering phase by det S E B , = exp(— 2iQ B ,(E)), where S E B , is the S-matrix with boundary 
condition B'. 
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where the A • B are the eigenvalues of — A n B , the Laplacian in fl with boundary conditions B. 
The S -matrix S E B , with boundary condition B' can also be described purely in terms of A z , 



Here 7 is the restriction to T, and 7^ is the normal derivative on T. Note that one of the 
consequences of Eq.(1.4) is a pointwise spectral duality, see Theorem 3.3 below. Finally, we 
describe an algorithm for computing the (^-function which is based on the double layer potential 



This paper can be read in two ways. On the one hand, it can be viewed as a collection of 
identities relating Dirichlet, and Neumann scattering to Dirichlet and Neumann membranes in 
terms of single and double layer potentials on the boundary of fi, through certain ^-functions. 
Many of these identities, although reminiscent of work by [KR, HS, SU, BS2] are in fact 
new. On the other hand, we give detailed proofs for those identities involving the double layer 
potential, and relating Dirichlet conditions on one side of the boundary to Neumann conditions 
on the other. 

Acknowledgments. Our interest in the problems discussed in this paper has been provoked by 
the papers of Smilansky et al. [SU] and Steiner et al. [BS2]. Part of our findings are a direct 
consequence of fruitful discussions with and inspiring seminars by members of these groups. 
In addition, we have profited from helpful discussions with V. Ivrii, A. Jensen, V.S. Buslaev 
and D.R. Yafaev. They all have contributed to clarify our views about the relevant issues. This 
work was supported by the Fonds National Suisse, and many of our contacts have been made 
possible by the semester "Chaos et quantification" at the Centre Emile Borel in Paris. 

2. Notations and Definitions 

We consider domains O which we call "standard domains." 

Definition. A domain il e R 2 is called a standard domain if 

a) il is bounded. 

b) T = <9fi is piecewise C 2 , with a finite number of pieces. 

c) The angles at the corners are non-degenerate, i.e., neither nor 2ir. 

d) The complement fl c of fl is connected. 



- It should be noted that the definition allows for domains which consist of several pieces. We 
shall, however only deal with the case of a connected domain to keep the notation simple. 

- The theory would be somewhat easier, with bounds which are not really any better, if we 
restricted our attention to smooth domains. However, in view of applications and examples, 
we think that the inclusion of corners is important. 




(1.4) 



Remarks. 



To formulate our results, we need to define the various spectral densities and scattering 
phase shifts. 
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Notation. We shall use throughout the subscripts D and N to denote Dirichlet and Neumann 
boundary conditions, respectively. We denote by B, B' a choice of boundary conditions among 
{D,N}. 

Definitions. We define here the quantities iV B , @ B ,. Let B, B' be boundary conditions in 
{D,N}. 

- The quantity N B (E) denotes the number of eigenvalues below E, counted with multiplicity, 
of — A n B . Here, A n B is the Laplacian in 1) with boundary condition B on T. 

- The quantity B , (E) is the total scattering phase for the scattering operator in O c , with 
boundary conditions B' on T. I.e., det S E B , = e ~ 2i6)B ' (s) . It is normalized to B ,(O) = 0, 
and it is defined as a continuous function of E. 

Our analysis will be based on a study of the single and double layer potentials which we 
define next. We denote by G the Green's function 

G, = —^— z , ;eC\R+. 



and the integral kernel which goes with it: 

G z (x,x') = l -H^\V~z\x-x'\) = l -J (V~z\x-x'\)--Y {V~z\x-x'\), z e C \ R+ . 

(2.1) 

These are, respectively, the Hankel and Bessel functions (in the notations of [AS]). Note that 
G z is the free Green's function, and the interaction will be described purely in terms of the 
boundary layer operators. Furthermore, the precise form of these functions is not relevant for 
our purpose, it suffices to know their asymptotic behavior for large and small arguments. 

Since we assume that the domain il is connected, we can parameterize the boundary by 
arclength, by a map s i-> x(s) mapping [0, lit) into R 2 . Here, we assume without loss of 
generality that the length of T is 2ir. 

We start by defining the "restrictions to the boundary." Let / be a function on R 2 . Then 

( 7± /)(a) = \imf(x(s)±eN(s)), 

(1nJ)(s) = ^N(s) ■ (Vf)(x(s)±eN(s)) . 

Here, N(s) denotes the outward normal to T at x(s) and ± indicates whether the limit is to be 
taken (along the normal) from the outside of il (+) or the inside (— ). Whenever the direction 
of the limit is irrelevant, we omit the index ±. In the corners, this definition is problematic, 
but since we only look at integral kernels, this does not matter. The following identities show 
where the various restrictions are defined: It is well known [Ne] that for all /? > -|, 

7 : < C (R 2 ) - L\T) , 

7* : £ 2 (H - ^(R 2 ) , 

7^ : ^(R 2 ) - £ 2 (r) , 

7=v : L\T) -ir^+^R 2 ), 
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from which the appropriate domains of 7 and 7^ can be read off. The notation is as follows: 
Let A = (1 — d 2 ) 1 ^ 2 , where d s is the derivative with respect to arclength on Y. Then, 

H p = {ueL 2 {Y) : A^u G L 2 (T)} , 

for P > 0, i?f omp is the subspace of functions with compact support, and H^ oc denotes functions 

which are locally in H 3 , see [H]. We also define i^(R 2 ) = {/ G L 2 (R 2 ) : (1 - A)^ 2 / G 
L 2 (R 2 )}. 

We next define the Single Layer Potential & 1 and the Double Layer Potential $ 2 . For 
x G R 2 \ F, we have 

= J ds 1 G z (x,x(s'))u(s') , 
($ 2 tt)(a:) = J ds'N(s') ■ {V 2 G z )(x,x(s))u{s') . 

The notation V 2 means the gradient with respect to the second variable of G z . (The z- 
dependence of ^- is implicit.) A more suggestive notation is 

*i = G z l* > $ 2 = ^ • 
Since G z maps # c "^ p (R 2 ) to H 2 ~/(R 2 ) for all 0, we see that 

<&! : L 2 (r) - < C (R 2 ) , $ 2 : L 2 (r) - ^(R 2 ) , 

for all S < §. 

We finally define 

B z = \{l+ + l-)(G zl * N ), (2.2) 
C z = l N {G zl * N ). 

Henceforth, we will omit the parentheses around G z . It is a well-known fact that the jump 
discontinuity of the single and double layer potential is 1: More precisely, we have the relation 

1 = (7+-7-)(Gj;) , 
-1 = (7jv + -7at_)(G,7*)- 

The operators A Z ,B Z , and C z are defined as maps between the following spaces: 



A„ : L 2 (T) -> ^(r) 

L 2 (r) -> L 2 (r) 
i? 1 ^) -> L 2 (r) 



z 

B 



z 

C, 
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For A z , this was shown in [EP1], for B z we will show it below, and for C z it will be a 
consequence of the bounds on A z and B z . 

The following expressions for the integral kernels may make explicit calculations more 
readable [CH]: 

A z (s,s') = ( 7 G zl *)(s,s') = G z (x(s),x(s')) , (2.4) 
B z (s,s') = ^N(s')-V 2 G z (x(s),x(s')), (2.5) 
Bj(s,s') = -^N(s)-V 1 G z (x(s),x(s')), (2.6) 

where V ■ is the gradient with respect to the j th variable. The operator B T is the transpose of 
B. The case of C is more complicated and will be handled in the Appendix. 

Finally, we define the main objects of this paper, namely (-functions, one for each of the 
boundary operators. 

Definition. We define 4 (-functions. We choose some negative number E and define, for 

zec\n + : 

Cdd(^) = det {A- E ]A Z ) , 
C DN (^) = &et({\-B Eo )-\\-B z )) , 
C ND (^) = det{{\ + B Eo )-\ l 2 +B z )) , 
C NN (^) = det {C E lC z ) . 

The first index will refer to the interior boundary condition and the second to the exterior 
boundary condition. 

Remark. There is a close relationship among the 4 (-functions, which is a consequence of the 
identity: 

C z = {\ + B T z )A-\\ - B z ) = (I - BF z )A-\\ + B z ) . (2.7) 

The identity Eq.(2.7) follows from the following considerations: Fix u and define the 
single layer potential $ = <^ x u. Then we have 7 $ = Au and 

lN _$ = (\ + B T )u. 

We can write the same function $ as a double layer potential: if $ = $ 2 w inside O, then 

7_$ = (\- B)v and 

1n _ $ = Cv . 

Hence we find 

Au=(±-B)v, (± + B T )u = Cv, 

from which the first identity in Eq.(2.7) follows. The second identity can be obtained by 
repeating the above arguments for fl c in place of ft, i.e., the limits are taken from the outside. 
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Using the identity Eq.(2.7), it is almost obvious that it suffices to study the ^-functions 
Cdd' Cdn' and Cnd- Then Cnn can be expressed as 



C NN (^) = det(C£C7j 

= det ((I - BzX'A^il + Slj- 1 ^ + BjK-^i - B z )) 

= det ((I + BlX\\ + SJJAJ^I - B z ){\ - B Eq )~ 1 A Eq ) 

= C ND (^) • det - - B^)" 1 ^) 

= CndOO • det {A E A-\\ - B z ){\ - S^)" 1 ) 

= CndOO • det {AzA- 1 ) • det ((I - B z ){\ - S^)" 1 ) 

= Cnd(^) • Cdd(^) • Cdn(^) • 



This identity allows to avoid the use of C which is more complicated to compute than A or B. 

3. The Relation Between the Scattering Phase and the Density of States 

Our main result is the following set of identities: 

Theorem 3.1. For any choice of B, B' £ {D, N}, the total scattering phases, the integrated 
density of states and the Q functions are related (for E > 0), by 



Remark. The relation Eq.(2.7) is reflected through the Eqs.(3.1), since the 4 possible left hand 
sides are linearly dependent. We furthermore have the bounds: 

Theorem 3.2. One has the following bounds, valid for E > 2: 



Remark. With slightly more complicated expressions due to threshold effects the formulas 
above extend to E = 0. 

Discussion. The above results describe a close relation between the integrated density of states 
and total scattering phase. Thus, they are much weaker than the spectral duality result ("inside- 
outside duality"), conjectured in [DS] and proved in [EP1], but they generalize and extend 
the pioneering result of [JK]. To complete the picture, we state here the result which relates 
individual eigenvalues and eigenphases: 



nN B (E) = e B ,(E)-ImlogC B) B'(£ + iO). 



(3.1) 



< 6 D (F) -vriV D (F) < const. E 1/2 log E , 
|0 N (F) — 7vN- D (E)\ < const. E 1/2 log E , 
|G D (E) -tt7V n (F;)| < const. E 1/2 log E , 

< tvN^(E) - 6 N (F) < const. E 1/2 log E . 



(3.2) 
(3.3) 
(3.4) 
(3.5) 
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Theorem 3.3. If fl is a standard domain, then E* is an eigenvalue of—A Q D of multiplicity m 
if and only m eigenphases of the S-matrix in Q, c with Dirichlet boundary conditions approach ir 
from below as E | E* . 

If fl is a standard domain, then E* is an eigenvalue of — A n N of multiplicity m if and only m 
eigenphases of the S-matrix in il c with Neumann boundary conditions approach from above 
as E { E*. 

Remark. The first part was shown in [EP1], the second part is new. We do not expect any 
similar result for the case when the boundary condition for the inside and the outside problem 
are not the same. Our convention of scattering phase is that the eigenphases of the (unitary) 
S-matrix are exp(-2i0 £ (.E)), I = 1,2, 

Remark. We next wish to comment on the bounds in Theorem 3 .2 and their possible optimality . 
The growth of Im log C, has, in our view, two different origins in the case of C DD an d Cnn 
when compared to the "mixed" cases £ DN and Cnd- ^ n me case of Cdd me growth can be 
traced back to the different Weyl asymptotics of D and iV D [SU], namely 

*N™(E) = ^E-^E^ + 0(1), 

e D (£) = Be + ^eV* + o(i), 

where the superscript W indicates the Weyl approximation, so that already at the "average" 
level the two quantities differ by {dfllE 1 ^ 2 . A similar formula holds for the pure Neumann 
case. On the other hand, in the case of Cnd' there is a different asymptotics since 

and therefore there is a cancellation of the terms of order E 1 ^ 2 at the Weyl level. 

We next discuss in detail the question whether this cancellation implies better bounds in 
Eqs.(3.3) and (3.4). (We neglect here the issue of eliminating the factor log E.) 

A first possibility might seem a proof using the properties of B. Indeed, in the integral 
kernel of B there appears the product N(s') ■ (x(s) — x(s')) /\x(s) — x(s') \ which goes to as 
s' — > 5, so that in the detailed bounds one additional order cancels when compared to the bound 
on d s A(s, s'), which occurs in the estimate for Im log £ DD (see Eq.(5.10)). However, we still 
cannot exclude that the bounds in Eqs.(3.3) and (3.4) are optimal, since actually the cancellation 
does not take place at intermediate distances, i.e., \s — s'\ = 0{1). 

A second possibility is provided by the very detailed results from the methods of pseudo- 
differential operators. The following discussion is a summary of the papers by Seeley, Melrose, 
Ivrii, Buslaev, Robert, and Vasil'ev, Safarov[VS]. The major new ingredient here is the notion of 
the set of periodic orbits of same length. We say that a billiard has property S (for synchronous) 
if there are "many" periodic orbits in the following sense. Consider a periodic orbit, of period 
T. Let (p T (z) denote the phase space point reached from z £ £1 x S 1 (initial position and initial 
direction of the orbit) after time T, i.e., the end point of the billiard trajectory (in phase space) 
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with initial time T. Let z* be the periodic point: ip T (z*) = z* . We say this orbit is "absolutely 
periodic" [VS] if 

f( Z ) EE W T (Z)-Z\ 

has a zero of infinite order at z = z* . In other words, the returning rays have infinite focusing 
in a neighborhood of z* . A billiard has property S if the set of absolutely periodic points has 
positive measure in phase space. (For example, if f(z) is a "devil's staircase," then it has a set of 
full measure of points where f(z) — f(z*) vanishes of infinite order when z — > z* .) Examples 
of billiards with property S are given in [VS]. If a billiard does not have property S, we say 
that it has property A (for asynchronous). When talking about scattering, this condition is to be 
applied to exterior orbits, which might for example be trapped in the "outside" of an obstacle. If 
a billiard has property S, the term of order E 1 ^ 2 in the Weyl series is modified by an oscillating 
amplitude, W(E), which can in principle be computed from the knowledge of the synchronous 
set. If fl is convex and the boundary is an analytic curve, then one has property A, but in most 
other cases, it is difficult to decide whether a domain has property A or S. The following table 
summarizes the known results for C°° boundary (in dimension 2, in odd dimensions, slightly 
more is known). The upper sign is for Dirichlet, the lower for Neumann boundary conditions. 



Property 


irN(E) 


Q(E) 


A 


\n\ ET \an\ E i/2 + ( E i/2 ) 


\n\ E± 1^1^/2 + ^1/2) 


S 


\n\ ET \^\ w(E)E i/2 + o{E i/2 ) 


l -^E + 0(E 1/2 ) 



If the boundary is Lipshitz then it is only known that Q(E) = \fl\E/4 + o(E), in all 
cases [Rl, R2]. Applying the results of this table to our questions, we see that the only known 
improvement over our bound seems to be a bound of o(E 1 ^ 2 ) in Eqs.(3.3) and (3.4) when 
property A holds. The disc is an example of this case [SU]. 

A third way to view these problems is in the context of scattering resonances [He] , although 
we have no rigorous results to offer. Consider domains with trapped orbits in VL C . Then we 
expect resonances and these may contribute to the growth of the ("-function. To decide how 
much they grow, one would have to know if these resonances stay near the real axis. If they do, 
they will contribute a term 0(E 1 ^ 2 ) to Im log Cdn- But ^ they move away fast enough from 
the real axis, as the energy increases, they might as well only contribute (9(1). In that case, the 
bound Eq.(3.3) would not be optimal. 
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4. The Krein Formula 

Using the methods of [EP2], one can derive from the identities of Theorem 3.1a corresponding 
set of Krein trace formulas and, in the case of identical boundary conditions only, a Green's 
formula. We shall state them here without proof. 

Definition. We let H = — A, and we define the "interacting" Hamiltonians 

^B,B' = ~~ ^fi,B © ~~ AfK,B' • 

Furthermore, we denote by A • B the eigenvalues (counted with multiplicity) of — A n B , and 
by S E B , the on-shell S-matrix for scattering on il c , with boundary conditions B'. With these 
notations, one has the 

Theorem 4.1. For every choice of B, B' e {D, N} and for all F £ 5(R + ) with support in 
{E : E > 0}, one has the identity 



The proof follows very closely the one in [EP2, Section 4], and will be omitted. 
We next state some identities for the S-matrix: 

Definition. We start by defining the operator H E , which is "restriction to the energy surface E." 
Let p e R 2 , p ■ p = E. Then 



Ti(F(H BtB ,) - F(H )) = E F (\b) 



n 



(4.1) 





R 2 
/•2-7T 



dPye-WrKv), 





E — 



Theorem 4.2. The S-matrix at energy E can be expressed as follows: 



1 



1 



1 



1 



2iriV E 'y*A E \ i0 'y?:* E 
27riS s7 *(I-i? s+ . )- 1 7 ^ + S* 
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Note that 7X^5 and all the other combinations of 7 and E used above are well-defined. 

We finally state two identities between Green's functions. Note that no such identity is 
available in the case when the inside boundary condition is not the same as the outside boundary 
condition. 

Theorem 4.3. The Green 's functions satisfy the following identities: 

G U u = G - G7M-SG , 
C NN = G - G7^C" 1 7atG ! • 

The boundary layer integral lends itself in a very natural and systematic way for the 
computation of CdN' as we ^ as f° r a determination of the eigenvalues and, to some extent of 
the scattering phases for the Dirichlet, resp. the Neumann problem. These algorithms work, at 
present, only for the case of domains where T is C 1 , i.e., corners are excluded, but jumps in 
the second derivative are allowed. Consider the integral kernel B(s, s'). To discretize it, we 
choose an ordered sequence of points s- on the boundary (this method is also used in [HS]). If 
the boundary is smooth, it is advisable to choose the points equidistant in arclength, since then 
the method is of infinite order in the step size [Ha]. In the other cases, we have chosen unequal 
steps, and in particular 2 points at distance at every discontinuity of the second derivative of 
the boundary, one point as the limit on either side. 

5. Bounds on the Double Layer Potential 

Our main ingredient for the proof of all the results stated so far are Structure Theorems, which 
describe the detailed regularity properties of the operators A z and B z , and hence, by Eq.(2.7), 
also those of C z . We first state this result: 

Definitions. If C is a compact operator, we let s n (C), n = 1,2,... denote the eigenvalues of 
(C*C) 1 ^ 2 in decreasing order. One defines the weak Schatten classes (for 1 < p < 00), as the 
set of those C for which 

{C) p = ^ V n^s n {C), 

n 

is finite. We also define the associated norms 

(00 
71=1 

We let P denote the orthogonal projection onto the constant functions in L 2 (Y). Recall also 
that A = (l-d 2 s ) 1/2 . 

Definitions. We need some precisions concerning the branch cuts in the definition of G z , 
cf. Eq.(2.1). Let £ denote {z : z e C, z £ R + }, and let 1Z denote the Riemann surface 
associated with the logarithm. The function is defined on 1Z, and the integral kernel G z is 
defined for z e £, with the convention that Imz 1 ^ 2 > for z e E. Then we have 
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Structure Theorem 5.1. 

- For z G 1Z the operator A A z has the following representation: 



AA Z = 1 + Q A + R A + T M, (5.1) 

where Q A is bounded and of norm \\Q A \\ < \, where R A is Hilbert-Schmidt, and where 
is trace class. Moreover, there is a constant K so that for z G S, one has the bounds 



(5.2) 



<Tj A > + iP log z) 2/3 < if|z| 3 / 4 |log*| , 
||Ti A )A+iP log^|| 2 < lf|*| 3 / 4 |log*|. 

The operators Q A and R A do not depend on z. 

For z ^TZthe operator B z has the following representation: 

B z = Q B + R B + jf ) , (5.3) 

where Q B is bounded and of norm \\Q B \\ < \, where R B is Hilbert-Schmidt, and where 
is trace class. Moreover, there is a constant K so that for z G £, one has the bounds 



(T { Z B % 3 < K\z\*'*\logz\, 
\\T^A\\ 2 < K|^| 3 / 4 |log^|. 



(5.4) 



The operators Q B and R B do not depend on z. 



Remarks. Note that B z already "contains" a derivative (the normal derivative), whereas in A z 
the derivative is provided by A. Note also that while the operator ^ seems to be absent from 
B z , it reappears naturally through the very definition of B z , Eqs.(2.2) and (2.3). So the results 
for A z and B z are in fact quite similar. In particular, as discussed before, we do not get better 
results for than for Tj a ^ , although such a result might have been expected from a local 
analysis. 



Proof. The proofs for the operators A z have been given in [EP2], [EP1], so we deal here only 
with B z . 

We note first that the Green's function equals G z (x, x ) = {i/A)H^\z 1/2 \x-x' |), where 
H { q ] is the Hankel function, cf. [AS, 9.1], 

H^\w) = J (w) + iY (w) . 

For z G £ (which is really a complex energy) we let k = z 1 ^ 2 . Starting from Eq.(2.5), we see 
that the integral kernel for B z is 

B z (s, s') = " . (5-5) 



Zeta Functions 

where G(x) = jH^\x). To study B z , we start with Eq.(5.5). First observe that 

d z Il£\z) = -(J^ + iY^z)) = -H[ L \z). 
Substituting the definition of G, we get 

B z (s,s) = -H[\k\x(s)-x(s)\) lx{s) _ x{sl)l ■ 
We introduce the notations 

B z (s,s') = l ±N{s')-D{s,s')H^\kr), 

with 

r(s,s') = \x(s) — x(s')\ , 
D(s,s') = (x(s) — x(s'))/r(s, s 1 ) . 



The local behavior of is given by 



H[ 1] (r) = -— + C(r(l + logr)) 



7rr 



and we introduce the regular part g of : 



g(r) = H[ 1 \r) + — 
7rr 

2i 

= — J\{f) logr + 0(r) 

7T 

= —r logr + 0(r) , 

7T 

cf. [AS, 9.1]. We will also need in the sequel the representation 

g '(r) = -logr + 0(1) = 9 -^- + 0{l). 
7r r 

With these notations, we define the regular and singular parts of B z . 

B^( s ,s') = ^N(s').D(s,s')g(kr(s,s')) 
so that B z = B sing + Bl eg . Note that B sing does not depend on z. 
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Proposition 5.2. The operator B z eg is in the Schatten class 2/3 and for all z e £ one has the 
bound 

(n /3 < const. |^| 3 / 4 (l + |log^|) . (5.7) 

Proof. In order to bound B z eg , we consider the quantity d s B z eg (s, s'). From Eq.(5.6), we get 
the representation 

g '(r) = -logr + C(l) = i^l + 0(l). 

7r ' r 

Furthermore, we have, with T(s) = d s x(s), 

d s r(s,s) = T(s)-D(s,s'), 

d s D(s,s') = -^—(T(s)-D(s,s')(T(s).D(s,s'))) , 
rys, s ) 

and therefore 

d s B r z eg (s,s') = l -^N(s')^(T - D(T ■ D))g(kr) + Dg\kr)krT ■ D^j 

= % -^N{s')- r {{T - D(T ■ D))g(kr) + D(^- + 0{\))krT ■ D^j 

= jiV(s') J (r(s)g(kr(s, s')) + krO(l)^j 

= k 2 (^N(s') ■ T(s) (± \og(kr(s, a')) + 0(l)j + 0(1)) 

= (t(s) ■ N(s') \og(kr(s, s')) + 0((kr) ) 

We next bound the Hilbert-Schmidt norm of d s B z eg and of B z eg . 

Lemma 5.3. There is a K < oo such that for all\z\ > 2, one has the bounds 

\\d s B z eg \\ 2 < if|z| 3 / 4 |logz|, (5.8) 
||^ eg || 2 < K\z\ 3/4 \\ogz\ . (5.9) 

Proof. We only show Eq.(5.8) and leave the proof of Eq.(5.9) to the reader. We split the integral 

|2 



J ds ds' \d s B z eg (s,s')\ 



into two parts, I = I < + I > corresponding to the integration region kr(s, s ) < 1 and its 
complement. We then get the bounds, observing that s and s' are bounded: 



J< = 0(k 4 ) f |l + |log(Ar 

Jkr<l 

< 0(k 4 )k~ 1 f dx(l + \logx\f < 0(k 3 ) 
Jo 



2 

(5.10) 
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We have used here that kr (s, s') < 1 implies | s — s' \ < 0(k~ 1 ). In the complement, we use that 
the integrand is bounded by k 4 (logkr) 2 , and a scaling argument shows that I > < 0(k 3 (logk) 2 ), 
as k — > oo. The proof of Eq.(5.8) is complete. 

Proof of Proposition 5.2. Note first that 

A = (l + (tfj 2 ) 1/2 = 11 + ^1 = U(l + d s ), 
where U is unitary. Therefore, the Schwarz inequality and Eqs.(5.8), (5.9), imply 

||A^ eg || 2 < K\z\ 3/4 \\ogz\ . 
Since the spectrum of A is + n 2 } neZ , we get (A -1 ^ < oo. We next use the inequality 

(C 1 C 2 ) p < 2 1 ^{C 1 ) q {C 2 ) r , 
valid for all q > 0, r > 0, p _1 = q~ l + r -1 , see [BS]. This implies 

(^ eg ) 3/2 = (A- 1 A J Br g ) 3/2 < 2 3 / 2 (A^) 2 (A~ 1 ) 1 < const. \\AB^\\ 2 . (5.11) 

The proof of Proposition 5.2 is complete. 

We consider next the operator B smg . Recall again that it does not depend on z. Our result 
here is 

Proposition 5.4. The operator B slng has the following representation: 

5 sing = Q B + R B , 
where Q B is bounded and of norm \\Q B \\ < \, and R B is Hilbert-Schmidt. 

Proof of the Structure Theorem 5.1. This is an immediate consequence of Proposition 5.4 and 
Proposition 5.2 and an inequality of the type of Eq.(5. 1 1). 

Remark. We shall use the notation to denote equality modulo Hilbert-Schmidt operators, 
i.e., s') « B 2 (s, s') means that B 1 — B 2 is Hilbert-Schmidt. 

Proof of Proposition 5.4. We shall localize the integral kernel of B smg near the corners and in 
the complement of this region, and give a different treatment for the two pieces. To localize, 
assume the boundary has exactly n corners, and define n smooth functions xf^ :Txr->R 
which are equal to 1 when s } s' are both close to the corner j, and such that they have disjoint 
supports. We then further decompose into the 2n functions 



opposite/ /\ (0)/ /\ opposite/ /\ 



V 7 J /VJ V 7 J /VJ,I 

xT me (s,s') = xf\s,s')-xtr(s,s') 
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where x° p poslte is 1 if s, s' are on opposite sides of (and close to) corner j, and otherwise, and 

^same = j _ ^Pposite _ Wg ^ describe ^ local behavior of B *™& Qn & sm0Qth part Q f p ^ 

near a corner of T. 

Lemma 5.5. On the smooth pieces of T one has 

B shlg (s,s') « -^-x(s) . (5.12) 

4"7T 

Here, x(s) is the curvature at s. 

Remark. These approximations hold near the corners as well if s, s' are on the same side of the 
corner. In particular, 

B™%s,s') h-^X? PP ° Site j (*>*') 



■±x(s) ^l-Exf P ° Site j (*,*') + O(s-s') 



From this and the compactness of T it follows that the corresponding contribution of 5 sing is 
Hilbert-Schmidt. 

Proof of Lemma 5.5. We recall the definition of B sing : 

Since we consider the smooth part of T, we get, with T(s) denoting the tangent vector, 

x(s) - x(s') = (s - s')T(s') - i(s - s') 2 x(s')N(s') + o((s - s') 2 ) , 

\x(s) - x(s')\ 2 = \s - s'\ 2 (l + o((s - s') 2 )) , 

\x(s) — x(s')\ = \s — s'\(l + o((s — s') 2 )) 5 

N(s') ■ (x(s) - x(s')) = - s') 2 x{s') + o((s - s') 2 ) . 

Substituting in Eq.(5.13), the assertion Lemma 5.5 follows. 

We next study the behavior "across" one corner, i.e., 5 sln g^PP° slte - 

Lemma 5.6. Modulo Hilbert-Schmidt operators, the operators B slng Xj PP ° slte are bounded in 

norm by -| | cos(^f)|, where a ■ is the interior angle at the corner j. Since the B slng Xj PP ° sl e have 
disjoint supports, it follows that there is a Hilbert-Schmidt operator K (which is independent of 
the energy z), for which one has the inequality 

n 

|| S sing y opposite _ R || < max ±| COs(^)| . (5.14) 
7— f j = l,...,n 

J = l 
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Proof. We assume that the boundary is locally straight. We leave to the reader the study of the 
Hilbert-Schmidt correction terms generated by curvature near a corner, see [EP1]. Assume for 
convenience that the corner is at the origin and use coordinates s and t = —s' near the corner, 
when s > and s' < 0. 




Fig. 1: The local coordinate system. 

Elementary geometry leads to 

\x(s) — x{s')\ 2 = s 2 + t 2 — 2stcosa , 
N(s') ■ (x(s) — x(s')) = N(s')x(s) = — ssino;. 

Redoing this calculation for s < 0, s' > 0, we find that for all s, s' satisfying ss' < one has 

N(s') • (x(s) — x(s')) = — |s|sina. 

We shall consider Q = B sing x opposite near one corner. We find, that Q equals (locally) 

1 



Q(s,s') = X (ss' <0) 



s sin a 



2tv s 2 + s' 2 + 2ss' cos a 



sign s 



47vi 



s + e %a s' s + e- %a s' 



We study, as in [EP1], the operator Q acting on L 2 (R). This is done by decomposing first 
L 2 (R) = L 2 (R + ) © L 2 (R + ), using the map u(s) i-> (u,(s),u_(s)) with 



u(s) 



u + (s), when 5 > 0, 
it_(— s), when s < 0. 



Having gone to unbounded coordinates, we can use them for an explicit calculation. 
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The operator Q is diagonalized by the Mellin transformation M, defined by 

1 f°° 

(Mf)(A) = -j=j^ dss^-^fis), 

as we shall show now. Indeed, this is intuitively clear since M. diagonalizes dilatations. Note 
that M : L 2 (R + ) — > L 2 (R) is unitary. With the above notation, we see that MQu is given by 



(MQu) a {\) 

= l r dss ix - i/2 r — (i - i ) u (s r ) 

where a e {+, — }. Replacing the integration variable s by ss' and noting that the integrand is 
homogeneous of degree iX + 1/2 in s', we get 

= -r^ a - i/2 (^-^i^)(^)(A), 

= c(A)(Mn_„)(A). 

Thus, Q becomes matrix multiplication under the Mellin transform. We next evaluate the 
integral c(A). Note that the integrand is 0(s~ 3 ^ 2 ) at infinity and 0{s~ 1 ^ 2 ) near 0. Therefore, 
for large R, we find 



The integrand is meromorphic in the annular sector {s : 1/R < \s\ < R, arg(s) e (0, 2ty)}. 
To evaluate the integral, we consider the contour given in Fig. 2. 

The integral over the circles which are concentric around the origin contributes 0(R~ 1 ^ 2 ) and 
the integral over the segment 1/R < s < R, arg(s) = 2tt equals (e 27 ")* A_1 / 2 c(A). Letting 
R — > oo, we obtain 

= -c(A) + (e 2 ™) iA - 1/2 c(A) 

JX-l/2 JX-l/2 

+ V Res(±- -). 

ze{e ia ,e i ( 2?r - a )} 

This leads to 

cosh((7r — a) A — ia/2) 
C(A) = 2cosh(7rA) ' 

Note that c(A) = when a = ir. 

We next observe that on the space L 2 (R + ) © L 2 (R + ), we have 
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Im s 



A 




Fig. 2: The contour used in evaluating the integral c(A) of Eq.(5.15). 

From this it follows at once that the spectrum of the corner contribution is 

= l cosh(A(7r- a) + jq/2) , A £ R} 
2 cosh(A7r) 

This set is shown in Fig. 3, for various values of a as a function of A. It is easy to check that the 
set cr(Q) is contained in the disk of radius \ cos(a/2). 

The "diagonal part" B sing is equal to B sing (l - ^ x JPP osite ). In this case, both s and s 1 
are on the same side of a corner, and we can reapply the bounds of Lemma 5.5. The proof of 
Proposition 5.4 is complete. 
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0.4 



0.2 




-0.2 



-0.4 

Fig. 3: The functions c(A) for angles a = 7rj/20, j = 1, . . . , 19. 

Appendix 

The definition and explicit study of the operator C z are more complicated than those of A z 
given in [EP1] and those of B z given in Section 5. These additional difficulties are generated by 
the appearance of two derivatives in the definition of C z , which make it an operator of order 1. 
However, when considering the integral kernel of C z , important cancellations of singularities 
take place, even in the presence of corners, and the main purpose of this appendix is a sketch of 
the mechanisms implying these calculations. Detailed notes on this problem can be requested 
from the authors. 

We fix first some simple notations. If x, y e R 2 , then we denote 

r = r(x,y) = \x - y\ , 

e = e(x,y) = V x r(x,y) = -V y r(x,y) = (x - y)/r(x, y) , 
x(s) = T(s) , 

f(s) = - K ( s )N(s) + Y / T j 6(s-s j ), 

j 

N(s) = x(s)T(s) + ^^6(8-8^. 

j 

The sums above are over the corners j, which are located at Sj G T, where 

T j =KmT(8 j + 6)-T(s j -6), 

and with a similar definition for iV". . The operator C z is defined starting from the double layer 
potential $ 2 m, which we defined for x £ T by 

($ 2 tt)(a;) = J ds N(s) ■ (V 2 G z )(x,x(s))u(s) . 
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Then 



(C z u)(s) = limiV(s) • V($ 2 u)(x(s) + eN(s)) 



for all s except the corners s = s •, j = 1, . . . n. It is not really obvious that the above limit 
exists, but we shall give a few arguments on the way to proving this. We want to argue only 
modulo Hilbert-Schmidt operators and we use the as sign to denote equality modulo Hilbert- 
Schmidt operators (with bounds independent of the energy). We note first the identity, with 

e = e(x, x(s')), r = r(x, x(s')), 

\7(® 2 u)(x) = z J ds' «(s')|e(e • N(s'))G(z 1/2 r) 

+ {2e(e.N(s'))-N(s'))^^ 

Here, again, G(x) = ^H^\x). In the limit process, the arguments of e and r will be 
(x(s) + eiV(s), x(s')), and then we find that (modulo Hilbert-Schmidt, as announced) 

f G(z 1/2 r) 
(C z u)(s) w lim / ds'u(s')N(s)-MJs,s')N(s')z V ' 
e|0 J T z L i z r 



-—lim / ds'u(s')N(s) ■ M(s,s')N(s')\ 



where M e is the 2x2 matrix 

2|e)(e|-l. 

Observing that the trace of M £ is 0, one can operate a certain number of cancellations. Further- 
more, it should be noted that the vectors T- and iV ■ describing the change of angle at a corner are 
orthogonal. This implies a cancellation of the first two orders of the most singular contributions 
at the corners, through the identity 

lim iV( 5j . -S)-M £ - T( Sj +6) = \im T( Sj - S) ■ M e ■ N( Sj + 6) . 

(All these statements hold modulo Hilbert-Schmidt contributions.) A lengthy calculation, 
involving the formula 

s^m = - (1+ ^.))iw-^w) 

shows, after integration by parts, that the result of these cancellations is 

(Cu)(s) w -—lim — — — - d I ds' u(s') e ' T ( S ^ 

V z A y 2tt £ |o 1 + ex(s) s 7 r V ' r 

w - lim d s J ds u'(s')G^z 1/2 r(x(s) + eiV(s), x(s')) 
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where u is the derivative of u. 

This result also allows us to compare C z to A z . In particular, it means that 

c z ~ i ~ 1 9 s A z i~ 1 d s , 
and therefore, by the Structure Theorem II of [EP2], we see that for all a E [0,1], 

(i - d 2 x a/2 c z {i - d 2 x /2 - 1/2 = \ + b + n + r z , 

where ||£?|| < |, 7i is Hilbert-Schmidt and T z is trace class. Note that B does not depend on 
z and that T z is actually in the Schatten class 2 /3. Of course, these facts also follow from the 
identity Eq.(2.7) and from the (explicitly proved) facts about A z , [EP2, Structure Theorem II] 
and B z , Structure Theorem 5.1. 
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